HEAT    TRANSFER   OF   FLOWING   GAS- SOLIDS 
MIXTURES   IN   A    VERTICAL   DUCT   AT 
DIFFERENT    TEMPERATURE    LEVELS 


by 


ERIK   WALDEMAR   STEPHANSEN,    JR. 

Diploma:     Dr.    Phil.    Almar  Ness  Technical  Institute 

Oslo,     Norway 


A   MASTER'S    THESIS 

submitted  in  partial  fulfillment  of  the 

requirements  for  the  degree 

MASTER  OF  SCIENCE 

Department  of  Mechanical  Engineering 


KANSAS   STATE    UNIVERSITY 
Manhattan,     Kansas 


1963 


Aoproved  by: 


Lv 

sn 

OoCUrne^i  TABLE   OF   CONTENTS 

INTRODUCTION 1 

Experimental  Work 1 

Theoretical  Work 3 

DESCRIPTION  OF  EQUIPMENT 3 

Heating  Section 4 

Adiabatic  Section 5 

Insulation 

Cooling  Section 5 

The  Cyclone 5 

Venturi 6 

Storage  Feed  Tank 6 

Feed  Mechanism 6 

Injector  Pump 7 

Entering  Section 7 

Recorder  7 

Power  Supply 7 

Meters 8 

CALIBRATION   OF   EQUIPMENT    23 

EXPERIMENTAL   PROCEDURE    33 

EQUATIONS 34 

RESULTS 46 

DISCUSSION   OF   RESULTS 48 

EXPERIMENTAL   ACCURACY 50 

RECOMMENDATION  FOR   FURTHER   WORK 50 

CONCLUSION    51 

ACKNOWLEDGMENTS 52 

NOMENCLATURE 5  3 

REFERENCES 54 

APPENDIX 55 


li 


INTRODUCTION 

The  study  of  heat  transfer  in  two  phase  flow  has  been  given  con- 
siderable attention  in  recent  years  because  of  the  promise  of  marked 
increase  in  the  heat  transfer  coefficient,    and  because  of  the  many 
fields  in  which  this  phenomenon  occurs.       The  nuclear  industry  is  very 
interested  in  increasing  the  heat  transfer  coefficient  of  gas-cooled 
reactors,    and  in  the  rocket  industry  (especially  that  dealing  with  solid- 
fuel  rockets)   complicated  two  phase  heat  transfer  situations  occur  in 
the  combustion  chamber  and  at  the  nozzle.       In  the  petroleum  industry 
this  phenomena  arises  in  some  cracking  processes,    and  it  is  believed 
that  this  industry  has  a  certain  amount  of  unpublished  literature  on  this. 

A  certain  amount  of  experimental  work  on  the  two  phase  flow  heat 
transfer  problem  has  been  published  to  date;    mostly  at  low  tempera- 
tures.     The  amount  of  theoretical  work  is  still  very  limited  because 
of  its  complexity. 

This  investigation  is  an  experimental  study  of  the  behavior  of  the 
local  heat  transfer  coefficient  for  different  air- solid  ratios  in  a  uni- 
formally  heated  tube  at  different  temperature  levels. 
Experimental  Work 

Systematic  experimental  work  was  first  reported  by  Farber  and 
Marley  at  the  University  of  California  in  1957.   (4)      They  passed  a  mix- 
ture of  air  and  aluminum  silica  catalyst  particles  (ranging  from  a  few 
microns  to  200  microns  in  size)  through  an  isothermal  heated  tube.     The 
average  Nusselt  number  was  reported  as  a  function  of  concentration  of 


particles  in  air.     For  mass  ratios  from  zero  to  one  the  Nusselt  number 
remained  almost  constant,    for  ratios  greater  than  one  the  Nusselt 
number  increased  steadily  and  had  doubled  itself  at  ratios  of  about  15, 
the  upper  limit  of  the  experiment.        This  work  was  later  extended  by 
Depew  (2)  who  studied  the  same  effects  by  using  uniformally  sized 
glass  beads  of  sizes  200,    70,    and  30  microns.       He  found  little  change 
in  Nusselt  number  as  a  function  of  concentration  for  the  bigger  parti- 
cles,   but  the  results  for  30  micron  particles  almost  matched  Farbar's 
earlier  curves.       Babcock  and  Wilcox  reported,    in  1959  ( l) ,    a  very 
extensive  study  of  a  two-phase  flow  experiment  carried  out  for  the 
nuclear  industry.      They  constructed  a  closed  loop,    in  which  a  mixture 
of  graphite  and  nitrogen  or  helium  was  pumped,    with  concentrations  as 
high  as   100  to  one.      They  obtained  an  increase  in  the  overall  heat  trans- 
fer coefficient  six  times  that  of  gas  alone.     In  i960,    Depew  (3)  reported 
on  further  work  in  this  field.     Special  attention  was  given  this  time  to 
the  change  in  local  Nusselt  number  as  a  function  of  concentration.     The 
study  was  carried  out  with  200  and  30  micron  glass  beads,   for  Reynolds 
numbers  13,  500  and  27,000.       For  the  200  micron  particles  he  found 
only  a  slight  increase  in  Nusselt  number  as  the  concentration  went  up 
while  for  the  fine  particles  (30  microns)  a  strange  decrease  was  observed 
in  the  Nusselt  number  to  about  20  per  cent  below  that  for  pure  air  for 
loadings  of  about  one,    followed  by  a  steady  increase  as  loading  increased. 
Quan  (6)  extended  this  work  by  using  different  particle  densities.     Lead 
and  glass  particles  of  ZOO  and  30  microns  were  used.     The  investigation 
showed  similar  results  as  Depew  for  glass  particles.       The  small  lead 


particles  behaved  similar  to  the  glass  particles,    but  the  decrease  in 
local  Nusselt  number  was  much  more  and  was  60  per  cent  below  that 
of  pure  air  at  its  minimum,    which  occurred  at  a  loading  ratio  of  1.2. 
Theoretical  Work 

The  only  known  theoretical  work  on  the  heat  transfer  of  a  flow- 
ing fluid-solid  mixture  was  published  by  Tien  in  l96l.    (l0)       He 
analyzed  the  problem  for  a  tube  with  constant  wall  temperature.      The 
analysis  is  very  interesting,    but  gets  extremely  involved  and  many 
simplifying  assumptions  were  made  to  obtain  a  solution.       The  results 
do  not  agree  well  with  Farbara's  and  Depwe's  experimental  results 
for  constant  wall  temperature.     Depew  (3)   extended  the  theoretical 
work  to  a  uniformally  heated  wall  and  found  qualitative  agreement 
with  his  experimental  work. 

DESCRIPTION  OF  EQUIPMENT 

The  equipment  was  constructed  in  such  a  way  that  the  particles 
would  pass  through  a  continuous  cycle  (closed  loop)  thereby  reducing 
the  number  of  particles  required.      The  loop  consisted  of  a  uniformally  % 
heated  heating  section  instrument  to  measure  the  local  heat  transfer 
coefficient  and  an  adiabatic  section  instrumented  to  measure  the  bulk 
temperature  and  to  study  the  mixing  length.       After  this  followed  a 
cooling  section  to  bring  the  mixture  to  approximately  room  tempera- 
ture before  the  particles  were  separated  from  the  air  in  a  high- efficiency 
cyclone.     At  the  air  outlet  of  the  cyclone  a  calibrated  venturi  was  in- 
stalled to  measure  the  air  flow.     The  particles  were  collected  in  a  delay 


tank  for  a  short  period  before  again  being  fed  into  the  loop.     The  air 
and  particles  were  pumped  around  the  loop  by  an  injection  pump 
driven  by  compressed  air.      The  pump  blew  straight  into  the  entrance 
section  to  create  a  good  mixture  of  air  and  particles  before  the  mix- 
ture entered  the  heating  section. 
Heating  Section 

The  heater  was  constructed  of  a  3/4-inch  diameter,    36  inch  long, 
stainless  steel  tube  with  an  inside  diameter  of  0.  710  inch.       The  wall 
thickness  was  24  gage.       For  recording  the  tube  wall  temperatures, 
30 -gage  thermocouples  were  spot  welded  to  the  tube  at  locations  shown 
in  Plate  II.       The   thermocouple  leads,    copper  and  constantan,    were  cut 
into  equal  lengths,    and  the  two  ends  were  gas  welded  together  into  a 
tiny  sphere.     The  sphere  with  the  leads  in  a  tangential  direction  to  the 
tube  were  then  carefully  spot  welded  to  the  tube.       Then,    electrical 
tape  was  placed  as  close  to  the  joint  as  possible  under  the  leads  and 
wrapped  around  the  tube  one  and  a  half  turns  covering  the  thermocouple. 
At  each  end  of  the  tube,    1  / 1 6  of  an  inch  from  the  end,    a  brass  flange, 
3/16  of  an  inch  thick,    was  silver  soldered  to  the  tube  for  connecting 
the  power  and  to  ensure  a  uniform  flow  of  current  into  the  pipe.     The 
leads  to  the  volt  meter  were  also  screwed  directly  onto  the  brass  flange 
to  insure  correct  voltage  readings  across  the  tube.       The  heating  section 
was  electrically  insulated  from  the  rest  of  the  system  by  two  reinforced 
fiber  bakelite  pieces  machined  as  shown  in  Plate  III.     Two  holes,    l/8-inch 
diameter,    were  drilled  perpendicular  through  each  block  to  provide 


pressure  taps  for  measuring  the  pressure  drop  across  the  heater. 
Adiabatic  Section 

The  adiabatic  section  was  stainless  steel  tubing,    36  inches  long. 
The  thermocouples  were  welded  to  the  tube  in  the  same  fashion  as  the 
heating  section  at  positions  shown  in  Plate  II. 
Insulation 

Both  the  heating  section  and  the  adiabatic  section  were  insulated 
in  the  same  manner.     Two  24-gage  galvanized  sheet  steel  tubes,    8 
inches  in  diameter  were  mounted  around  each  section.     When  every  part 
of  the  loop  had  been  tested  out,    the  containers  were  filled  with  vermi- 
culite. 
Cooling  Section 


The  entire  length  of  pipe  between  adiabatic  section  and  cyclone 
was  made  into  a  cooling  section.      Three  pieces  of  copper  pipe,    3/4  of 
an  inch  in  diameter  with  water  jackets,    were  used.     When  constructing 
the  equipment,    it  was  not  known  how  long  the  copper  pipe  bends  would 
last  against  the  eroding  effect  of  the  particles.       An  inspection  of  the 
pipes  after  completion  of  the  experiment  indicated  little  wear. 
The  Cyclone 

The  cyclone  was  constructed  in  two  parts,    the  upper  cylindrical 
housing  was  made  out  of  steel  to  prevent  wear.     The  conic  bottom  part 
was  made  out  of  glass  for  visible  inspection.     The  two  parts  were  glued 
together  with  Epoxy  resin.     The  cyclone  proved  to  be  almost  100  per 
cent  efficient  for  30  micron  particles.     For  details,    see  Plate  IV. 


Venturi 

The  venturi  was  made  of  3/4-inch  diameter  copper  pipe  by- 
forming  it  to  the  proper  dimensions  on  a  lathe.       The  orifice  had 
an  area  one  half  the  tube  diameter  which  proved  sufficient  for  ac- 
curate readings  on  the  manometer.       The  venturi  was  calibrated  by- 
passing a  steady  flow  of  air  through  it  and  letting  the  air  fill  a  gas 
meter  calibration  tank.       The  air  temperature  and  the  time  for  fill- 
ing the  tank  were  measured  together  with  the  pressure  drop  across 
the  venturi.     (See  calibration  curve,    Plate  VI.) 
Storage  Feed  Tank 

The  storage  feed  tank  was  constructed  out  of  a  poly  vinyl  flex- 
ible tubing     1-1/2  inches  in  diameter,    and  1  foot,    6  inches  long.       A 
flexible  connection  was  installed  between  the  cyclone  and  the  feed 
nozzle.       This  was  vibrated  to  insure  a  uniform  flow  of  particles. 
Unfortunately,    this  principle  gave  unsteady  flow  conditions  and  was 
changed  to  a  motor-driven  feed  mechanism  which  proved  more  sat- 
isfactory. 
Feed  Mechanism 

A  3/4-inch  diameter,    2-inch  long  copper  tube  was  split  3/  1 6  inch 
wide  and  I  inch  deep  at  one  side,    and  the  tube  was  placed  in  the  bottom 
of  the  feed  tank.     A  rubber  ball,    0.68  inches  in  diameter,   was  placed 
inside  the  tube  mounted  on  a  3/l6-inch  shaft.       The  ball  had  18  blades 
cut  axially  out  of  it.       The  shaft  was  connected  to  an  electric  motor, 
whose  speed  was  controlled  by  a  variac. 


Injection  Pump 

The  injection  pump  was  constructed  from  a  3/4-inch  diameter 
copper  elbow  through  which  a  I /4-inch  copper  pipe  air  jet  nozzle  was 
mounted.      The  jet  blew  into  a  venturi  also  formed  from  a  3/4-inch 
diameter  copper  pipe.       The  vacuum  created  on  the  upstream  side 
proved  adequate  to  transport  particles  from  the  feed  inlet  to  the  pump. 
Entering  Section 

The  entering  section  was  made  of  the  same  stailness  steel  pipe 
as  the  heating  section  (35  inches  long).       Three  inches  from  the  top, 
a  thermocouple  was  soldered  to  the  pipe  to  measure  the  entering 
mixtures  temperature.       The  thermocouple  joint  was  insulated  from 
the  surroundings  with  one-inch  glass  wool  insulation  covering  6  inches 
of  the  tube. 
Recorder 

A  Honeywell  high  speed  electronic  recorder  with  variable  sensi- 
tivity and  range  from  one  to  21  millivolts  full  scale  was  used.       The 
recorder  was  calibrated  for  each  set  of  runs  using  a  portable  millivolt 
potentiometer.       For  converting  the  recorder  readings  to  temperatures, 
the  National  Bureau  of  Standards  Reference  Table  for  Thermocouples, 
Circular  56l,    was  used. 
Power  Supply 

An  electric  welder  was  used  for  energy  control  to  the  heating  sec- 
tion.    For  the  highest  temperature  levels,    two  welders  were  used  in 
parallel. 


Meters 

A  General  Electric  precision  volt  meter  with  6-inch  full  scale 
deflection  was  used.       The  meter  had  two  scales,    0-3  volts  and 
0-15  volts.       Resistance  across  the  meter  was  650ft     for  the  0-3  volts 

scale. 

For  current  measurement,    a  50  millivolt  full  scale  volt  meter 
was  used  across  a  300  ampere  shunt  calibrated  for  50  millivolts. 
Full  scale  deflection  of  the  meter  was  6  inches. 

A  Meriam  manometer  with  one-hundreths  of  an  inch  divisions 
was  used  to  measure  the  pressure  drop  across  the  venturi. 


EXPLANATION   OF   PLATE   I 

Schematic  diagram  of  the  experimental  equipment 
to  study  heat  transfer  phenomenon  of  fluid  solid- 
mixtures. 
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EXPLANATION  OF   PLATE   II 
Thermocouple  positions  for  the  heating  and  adiabatic  section. 
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EXPLANATION  OF   PLATE  III 

Detailed  drawing  of  the  connection  between 
entering  and  heating  sections. 


PLATE      III 
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EXPLANATION   OF   PLATE   IV 


Detailed  drawing  of  the  cyclone. 
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EXPLANATION   OF   PLATE    V 


General  arrangement  of  apparatus. 


EXPLANATION   OF   PLATE    VI 

View  of  heating  section  with  instrumen- 
tation before  being  insulated. 
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EXPLANATION   OF   PLATE    VII 

Close  view  of  cyclone,    delay  tank,    feed  mechanism, 
injector  pump  and  entering  section. 
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CALIBRATION   OF   EQUIPMENT 

The  heat  loss  from  the  heating  section  was  measured  by  taking 
off  the  entrance  section  and  closing  the  bottom  of  the  heater.     A  low 
current  was  then  passed  through  the  section  for  about  4  hours  to  in- 
Sure  equilibrium  with  the  surroundings.     The  temperatures  of  the 
tube  and  the  heat  flux  were  recorded.      This  was  carried  out  a  four 
different  temperature  levels.     Total  heat  loss  was  then  plotted  against 
the  temperature  of  thermocouple  11,    see  Plate  V.       The  relationship 
gave  a  straight  line  as  expected.       It  was  decided  that  thermocouple  11 
would  give  the  best  reference  point  for  the  heat  loss  from  this  section. 
The  voltage  increase  across  the  tube  for  a  constant  current  was 
checked  by  using  the  thermocouples  as  contact  points.     It  was  assumed 
that  this  increase  should  be  a  linear  function  of  tube  length.      The  devia- 
tion is  shown  on  Plate  VII,    Fig.    I,    as  a  percentage  of  total  voltage  drop 
across  the  tube.       Notice  the  fairly  constant  deviation  of  0.  5  per  cent 
between  thermocouple  7  and  16.       About  half  of  this  can  be  accounted 
for  as  tube  expansion.       The  rest  might  be  caused  by  inaccurate  read- 
ings of  meter  or  a  change  in  tube  wall  thickness.       The  irregularity  of 
junction  14  is  believed  to  be  due  to  a  decrease  in  tube  wall  thickness. 

The  heating  section  was  then  tested  with  pure  air  at  the  275°  tem- 
perature level  to  check  the  local  Nusselt  number  at  x/D  =  46.8  against 
Sparrow's  (7)  analytical  expression. 

0.  77 
Nu      =     0.0245  (Re)  Pr  =  0.7 

The  resulting  data  checked  out  to  within  1.  5  per  cent  of  this  equation. 
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The  current  passing  through  the  heating  section  proved  to  affect 
some  of  the  thermocouples.       It  was,    therefore,    decided  to  check  this 
effect  for  each  temperature  level.       The  method  used  was  to  flip  off 
the  current  for  a  short  period  and  record  the  change  in  temperature. 
On  an  average,   the  change,   if  any,  was  within  1-2  °F,   except  for 
thermocouple  8  which  had  a    4    degree  depression  of  the  temperature. 
The  recorder  response  to  full  scale  deflection  was  one  second. 

To  study  the  grinding  effect  on  the  particles,    a  sample  of 
300  grams  was  taken  and  passed  through  the  loop  for  one  hour  at  a 
feed  rate  of  unit  mass  ratio  of  air  and  solids.       It  was  found  that 
l/3  of  the  particles  had  received  small  fractures  on  their  surface, 
and  a  few  out  of  a  thousand  had  cracked. 


EXPLANATION  OF  PLATE    VIII 

Heat  loss  from  heating  section  as  a  function 
of  temperature  recorded  by  thermocouple  11. 
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EXPLANATION   OF   PLATE   IX 
Calibration  curve  for  venturi. 
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EXPLANATION   OF   PLATE   X 

Fig.    1.     Percentage  deviation  of  the  measured  voltage  drop 
across  tube  compared  with  an  assumed  linear  drop. 

Fig.    2.     Temperature  profile  across  tube  without  air  flow. 
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EXPLANATION  OF   PLATE   XI 
Size  distribution  curve  of  the  30  micron  glass  beads. 
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EXPERIMENTAL   PROCEDURE 

The  equipment  was  always  started  up  at  least  two  hours  prior 
to  any  run  to  insure  equilibrium  conditions  for  all  parts.       This 
period  was  used  to  adjust  air  flow,    compressor  bleed-off,    and  heater 
Current  for  the  desired  temperature  level.       The  cold  junction  was 
filled  with  ice  and  the  recorder  calibrated.       Thermocouple  16, 
x/D  -  46.  8  diameters,    was  used  as  temperature  level  reference  point 
and  kept  within  a  few  degrees  of  the  selected  level. 

When  the  air  flow  and  all  temperatures  had  stayed  constant  for 
about  half  an  hour,    the  first  set  of  readings  were  taken  for  pure  air. 
Particles  were  then  fed  into  the  apparatus  and  pumped  around  the 
loop  to  the  delay  tank.       The  loop  was  operated  for  an  additional  ten 
minutes  to  allow  conditions  to  stabilize.       The  particles  were  then  fed 
into  the  loop  and  the  air-flow  and  heat  flux  adjusted.       This  took  from 
two  to  five  minutes.       At  high  particle  concentrations,    it  was  found 
harder  to  adjust  the  system.       Another   15  minutes  was  allowed  to  reach 
equilibrium.       The  temperature  of  thermocouple  16  was  then  recorded 
followed  by  readings  of  the  other  thermocouples  from  1  to  25,    and  then 
back  to  16  to  check  for  drift,    since  this  took  at  least  five  minutes.       If 
everything  looked  fine,    the  current,    voltage  and  pressure  drop  across 
the  heater  were  recorded  and  the  feed  rate  was  measured.       To  measure 
the  feed  rate,    a  beaker  was  passed  under  the  feed  mechanism  and  the 
sampling  time  recorded  with  a  0.  2  second  division  stop  watch.       At  low 
feed  rated,    below  unit  mass  ratio,    the  same  sample  was  used  a  few 
times,    but  at  high  feed  rates,    a  new  sample  of  particles  was  used  for 
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each  run  to  insure  no  change  in  particle  quality. 

EQUATIONS 

The  main  purpose  of  this  investigation  was  to  study  the  change 
in  Nusselt  number  as  a  function  of  mass  flow  ratio  for  different  tem- 
perature levels.       Three  different  levels  were  chosen:     185,    275,    and 
390°F. 

The  equations  and  relationships  used  are  listed  below. 

The  local  Nusselt  number  is  taken  as 

h     D 

Nu      =     — ^ 

x  k 

where  the  local  heat  transfer  coefficient  is 


h*    °    T        -  V 

w  bm 

x  x 

and  the  average  mixture  temperature  is  given  as 

T  T    +  ttD    qx 

bm  ll  +       WC+   W     C 

X  s      s 

All  properties  were  evaluated  at  the  local  bulk  temperature  except 
for  specific  heats,     C  and  C    ,   which  were  evaluated  at  I00°F.   (Change 
of  specific  heat  was  less  than  one  per  cent  for  temperature  range  used.) 

An  artificial  Reynolds  number  was  used,    assuming  the  mixture 
properties  were  the  same  as  for  pure  air.       Total  mass  air-flow  was 
kept  constant  throughout  the  whole  experiment.       This  gave  considerable 
variation  of  Reynolds  numbers  for  the  different  temperature  levels. 
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The  local  Reynolds  number 

P  Du 
Re  =    


or 

4W 


Re  = 


tt[i.  f) 


The  mass  of  air  and  particles  were  given  as  a  mass  ratio, 

W 


W 
and  axial  distance  measured  from  air  inlet  side  of  heater  as 


x 

T7 


EXPLANATION   OF   PLATE    XII 

Local  Nusselt  number  variation  with  mass  ratio  for  the  three 
different  temperature  levels.     Mass  flow  of  air  constant  at 
29  lb.    per  hr. 


PLATE     XII 
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EXPLANATION   OF   PLATE   XIII 


A  general  picture  of  the  variation  of  the  Local  Nusselt  number 

in  relation  to  recorded  and  computed  values,    all  as  a  function 

of  mass  ratios.     Temperature  level    -     185    F. 

Fig.    1.     Recorded  bulk  temperature  and  computed 
bulk  temperature. 

Fig.    2.     Local  Nusselt  number  at    x/D  =  46.8. 

Fig.    3.     Net  heat  flux  to  tube. 
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EXPLANATION   OF   PLATE   XIV 


A  general  picture  of  the  variation  of  the  local  Nusselt  number 

in  relation  to  recorded  and  computed  values,    all  as  a  function 

of  mass  ratios.      Temperature  level    -     275     F. 

Fig.    1.     Recorded  bulk  temperature  and  computed 
bulk  temperature. 

Fig.    2.     Local  Nusselt  number  at    x/D  =  46.8. 

Fig.    3.     Net  heat  flux  to  tube. 
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EXPLANATION   OF   PLATE   XV 

A  general  picture  of  the  variation  of  the  local  Nusselt  number 

in  relation  to  recorded  and  computed  values,     all  as  a  function 

of  mass  ratios.     Temperature  level    -     380     F. 

Fig.    1.     Recorded  bulk  temperature  and  computed 
bulk  temperature. 

Fig.    2.     Local  Nusselt  number  at    x/D    =    46.8. 

Fig.    3.     Net  heat  flux  to  tube. 
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EXPLANATION   OF   PLATE   XVI 

A  comparison  of  the  present  experimental  work  with 
other  experimental  and  theoretical  work. 
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RESULTS 

The  main  purpose  of  this  experiment  was  to  determine  the 
variation  of  the  local  Nusselt  number  as  a  function  of  particle  load- 
ing for  different  temperatures  in  the  heating  section. 

The  results  of  the  investigation  have  been  plotted  against  dif- 
ferent parameters  to  show  effects  of  interest  and  to  compare  them 
with  other  investigations. 

In  Plate  XII,    the  results  for  the  three  temperature  levels  are 
shown.       One  notices  immediately  the  great  decrease  in  Nusselt 
number  of  about  20  -   25  per  cent  as  the  mass  ratio  increases  from 
zero  to  0.7,    followed  by  a  rapid  increase  until  a  mass  ratio  of  about 
two  is  reached.       From  here,    it  seems  to  level  out  or  even  decrease 
as  the  mass  ratio  approaches  4.0.       Four  was  the  upper  limit  of  the 
feed  mechanism. 

The  Nusselt  number  shows  little  change  between  the  three 
temperature  levels  as  a  function  of  particle  concentration.       The  change 
in  the  total  value  of  the  Nusselt  number  is  caused  by  a  change  in 
Reynolds  number  ,    due  to  changes  in  the  mean  terminal  temperature 
of  the  mixture. 

In  Plates  XIII,    XIV,    and  XV  the  local  Nusselt  number  is  shown 
together  with  the  adiabetic  section's  bulk  temperature  and  the  heat 
flux  for  different  load  ratios.       The  outlet  temperature  of  the  heating 
section  is  believed  to  represent  the  bulk  outlet  temperature  of  the 
heating    section  to  within  a  few  degrees.     Notice  the  steady  drop  in 
bulk  temperature  as  the  loading  ratio  goes  to  0.7,    followed  by  a  slow 
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but  constant  temperature  drop  all  the  way  to  a  loading  ratio  of  four. 
The  calculated  bulk  temperature  given  much  the  same  shape,    and 
almost  corresponds  at  high  solid  loading  ratios.       This  is  a  good 
indication  of  the  accuracy  of  the  measured  feed  rate  of  particles. 

For  pure  air,   the  computed  bulk  temperature  i§  almost  10    F 
higher  than  the  measured  one  which  indicates  that  the  measured  air 
flow  might  be  too  low.       An  increase  of  3.  5  per  cent  of  the  air  flow 
would  correct  this  difference.       Many  tests  were  made  to  verify  this, 
but  without  success. 

In  studying  the  various  results  shown  in  Plates  XI,    XII,     and 
XIII,  it  appears  as  if  there  are  three  identifiable  regions.       In  the 
first  part,    almost  no  change  in  heat  flux  takes  place  but  a  disting- 
uishable drop  in  mixture  bulk  temperature  is  present.       In  the 
second  region,    the  heat  flux  increases  very  rapidly.       Finally,    a 
region  in  which  both  the  heat  flux  and  bulk  temperature  remain 
fairly  constant.       The  reason  for  these  three  regions  is  believed  to 
be  due  to  particle  concentration  and  its  interference  with  the  bound- 
ary layer. 

The  results  for  the  275     F  temperature  level  has,    on  Plate  XV, 
been  compared  with  other  investigator's  results.       A  general  agree- 
ment is   seen  with  both  Depew's  (5)   and  Quan's  (4)   experimental  values 
for  mass  ratios  up  to  two.       Between  1.5  and  3,    most  experimental 
values  seem  to  agree  closely  with  Depew's  (5)   theoretical  values. 
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DISCUSSION   OF   RESULTS 

It  might  be  worth  noting  the  good  agreement  for  the  local 
Nusselt  number  at    x/D  =  46.8  for  pure  air  with    Sparrow's  (7) 
analytical  equation  which  gave  an  agreement  of  5  per  cent  for  the 
380     F  temperature  level,    and  1.5  per  cent  for  the  275°  F  level. 
For  the   185     F  temperature  level  the  agreement  was   15  per  cent. 

In  Plate  XIV,   where  the  experimental  results  for  the  275°  F 
temperature  level  is  compared  with  other  investigations,    it  is  in- 
teresting to  note  the  general  agreement  of  the  shape  of  the  curves, 
although  the  actual  values  differ  considerably.       Quan's  results 
are  based  on  a  higher  Reynolds  number  (l5,  000),   which  accounts 
for  the  higher  values  of  Nusselt  number    at  lower  concentrations. 
The  experimental  results  for  all  three  investigations  are  fairly  close 
between  mass  ratio  1.  5  and  3,    and  agree  also  with  Depew's  theoret- 
ical curve.       It  should  be  mentioned  that  Depew's  theoretical  func- 
tion predicts  a  greater  drop  in  the  local  Nusselt  number,    at    mass 
ratios  of  one,    for  large  particles  than  for  small  particles,   while 
experimental  work  (3),   (6)  shows  the  opposite  trend. 

In  studying  Plates  XIII,  XIV,  and  XV,  the  following  conclusion 
was  drawn  to  explain  the  local  Nusselt  number  behavior  for  different 
parucle  to  air  concentrations. 

The  definition  of  the  heat  transfer  coefficient  between  tube  and 
fluid  is  defined  as  the  heat  flux  of  the  tube  divided  by  the  difference 
between  tube  wall  temperature  and  the  fluid  bulk  temperature.     But 
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the  main  resistance  for  heat  flow  between  tube  and  fluid  is  the  tempera- 
ture gradient  at  the  wall,    in  the  laminar  sub-layer.       If  a  few  particles 
are  added  to  the  fluid,    their  interference  with  the  boundary  layer  will 
be  very  limited.       But  they  will  cause  a  definite  change  in  the  mix- 
ture bulk  temperature. 

While  the  heat  flux  is  expected  to  stay  almost  constant  because 
of  the  small  interference  with  the  boundary  layer,    both  of  these 
tendencies  are  clearly  seen  from  the  recorded  heat  flux  and  bulk 
temperature  for  loadings  up  to  0.7.       As  the  number  of  particles  in- 
creases,   they  will  start  to  cause  an  appreciable  disturbance  of  the 
boundary  layer,    and  increase  the  temperature  gradient  at  the  wall. 
If  the  tube  wall  has  to  be  maintained  at  a  constant  temperature,    the 
heat  flux  must  be  increased  to  compensate  for  the  increased  heat 
transfer.       This  tendency  seems  definitely  to  be  the  case  for  the  heat 
flux  of  the  tube  for  loading  ratios  from  0.  7  to  2.  5.       If  this  disturb- 
ance of  the  boundary  layer  is  a  continuous  proves s  for  higher  concen- 
trations,   it  is  hard  to  determine  from  this  experiment.     More  data 
are  needed  to  predict  what  mechanisms  might  be  involved  for  higher 
concentrations . 
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EXPERIMENTAL   ACCURACY 


The  accuracy  for  each  measured  quantity  has  been  listed  below. 


Quantity 

Base 

Value 

Estimated  Deviation, 

36 

per  cent 

L    inch 

0.1 

0     inch 

0.  710 

1.0 

x    inch 

20 

0.  2 

Voltage    volt 

3.00 

0.5 

Current    amp 

100.00 

0.  5 

T,      °F 

80 

1.0 

T         °F 
w 

W      Ib/hr 

300 
29 

1.0 
1.0 

Feeding  time, 

sec. 

10 

2.0 

Mass  of  particles, 

grams  100  0.01 

The  estimated  accuracy  for  the  local  Nusselt  number  should  be  7 
per  cent. 

RECOMMENDATIONS   FOR  FURTHER   WORK 

1.  It  would  be  of  great  interest  to  see  if  the  local  Nusselt  number 
would  behave  the  same  for  particles  smaller  than  30  microns.     All 
other  parameters  should  be  maintained  as  for  this  work. 

2.  To  study,    experimentally,    the  particle  distribution  and  the 
velocity  profile  in  a  vertical  pipe  with  turbulent  flow.       This  is  a  very 
complex  study,    and  for  this  reason  few  investigators  have  tried  it. 
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CONCLUSIONS 

The  statements  made  pertain  to  a  uniformly  heated  tube  through 
which  air  and  particles  are  passed  having  a  fully  developed  velocity 
profile. 

1.  The  temperature  of  the  tube  wall  has  little 

effect  on  the  local  Nusselt  number  for  the 

three  temperature  levels  tested,    (l85    F,    275    F 

and  380°  F  ). 

2.  The  decrease  in  local  Nusselt  number  followed 

by  an  increase  as  concentration  of  particles 
increases,    can  be  shown  to  be  related  to  the 
definition  of  the  heat  transfer  coefficient  and 
to  the  particles  increasing  interference  with 
the  boundary  layer  as  concentration  goes  up. 

3.  The  local  Nusselt  number  as  a  function  of  particle 

concentration  agrees  fairly  well  with  the  results 
of  other  investigators. 
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NOMENCLATURE 

C  =  Heat  capacity  of  air  BTU/lb  F 

C  =  Heat  capacity  of  particles     BTU/lb  F 

D  ■  Tube  diameter,  ft 

h  =  Tube  wall  heat  transfer  coefficient  BTU/hr  sq  ft  F 

k  =  Air  therminal  conductivity    BTU/hr    ft  F 

N  =  Nusselt  number 
u 

q  -  Tube  wall  heat  flux    BTU/hr,    sq  ft 

R  =  Reynolds  number 

Ti  =  Mixture  entering  temperature,      F 

T,  =  Mixture  bulk  temperature,      F 

T  =  Tube  wall  temperature  ,      F 

w  r  ' 

W  =  Mass  flow  of  air,   Ib/hr 

W  -  Mass  flow  of  particles,     lb/hr 

P  =  Gas  density,     lb/ft3 

p.  =  Viscosity  of  air,     lb/ft  sec. 

x  =  Denotes  local  position  along  tube 
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In  recent  years,    considerable  attention  has  been  given  to  the 
study  of  two-phase  flow  with  heat  transfer  resulting  from  interest 
in  many  fields  such  as  the  nuclear  industry  and  rocket  industry. 
The  amount  of  thorough  experimental  and  theoretical  work  is  still 
fairly  limited  and  much  more  is  required  to  get  a  better  under- 
standing of  the  phenomena. 

This  investigation  is  an  experimental  study  of  the  problem 
carried  out  in  a  uniformly  heated  vertical  tube.      The  apparatus  was 
constructed  in  the  form  of  a  loop  through  which  a  mixture  of  air  and 
glass  beads  were  pumped.     It  consisted  of  an  entering  section  to 
obtain  a  fully  developed  velocity  profile,    a  heating  section  instru- 
mented to  measure  the  local  heat  transfer  coefficient,    an  adiabatic 
section  to  measure  the  mixture's  bulk  temperature  and  mixing  length, 
a  cooling  section  and  a  cyclone  to  separate  air  and  particles  in  order 
to  measure  these  quantities.     An  air  injection  pump  was  used  to  pump 
the  mixture  around  the  loop.     The  air  flow  was  kept  constant  at  29  lb/hr 
throughout  the  experiment.       The  heating  section  was  operated  at  three 
temperature  levels  -   185,    275,    and  380°  F,    and  the  particles  were 
spherical  glass  beads  30  microns  in  diameter. 

The  results  showed  that  the  local  Nusselt  number  behaved  much 
the  same  at  all  three  temperature  levels,    but  behaved  rather  peculiar- 
ly as  the  number  of  particles  increased.       The  local  Nusselt  number 
first  decreased  as  the  air  to  particle  mass  ratio  increased,    and  at 
mass  ratio  0.  7  it  was  20-25  per  cent  below  that  for  pure  air.     Nusselt 


number  then  increased  again  rapidly  as  mass  ratio  was  further 
increased,    but  at  a  ratio  around  two  started  to  flatten  out  and 
even  decreased  for  ratios  close  to  four.       The  reason  for  these 
three  distinguishing  regions  is  believed  to  be  due  to  three  dif- 
ferent particle  motions  and  their  interference  with  the  boundary 
layer  and  velocity  profile. 

The  results  for  the  275°  F  temperature  level  were  compared 
with  the  results  of  other  investigators,    and  a  general  agreement 
was  found  for  the  local  Nusselt  numbers  for  mass  ratios  up  to  two. 
For  higher  loadings,    the  experimental  results  agree  somewhat 
better  with  the  theoretical  values  obtained  by  Depew. 


